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We discuss the features of nonequilibrium growth problems, their scaling description and their 
differences from equilibrium problems. The emphasis is on the Kardar-Parisi-Zhang equation and 
the renormalization group point of view. Some of the recent developments along these lines are 
mentioned. 



I. INTRODUCTION 

How to characterize the degree of roughness of a sur- 
face as it grows and how the roughness varies in time have 
evolved into an important topic due to diverse interest in 
physics, biology, chemistry and in technological applica- 
tions. One crucial aspect of these nonequilibrium growth 
processes is the scale invariance of surface fluctuations 
similar to the scale invariance observed in equilibrium 
critical point phenomena. Although different kinds of 
growths may be governed by distinct natural processes, 
they share a common feature that the surface, crudely 
speaking, looks similar under any magnification and at 
various times. This nonequilibrium generalization of scal- 
ing involving space and time (called "dynamic scaling" ) 
makes this subject of growth problems important in sta- 
tistical mechanics. 

Growth problems are both of near-equilibrium and 
nonequilibrium varieties, and, therefore, they provide us 
with a fertile ground to study the differences and the ex- 
tra features that might emerge in a nonequilibrium situ- 
ation 1^,1^ . Take for example the case of crystal growth. 
In equilibrium, entropic contributions generally lead to 
a rough or fluctuating surface, an effect called thermal 
roughening, but, for crystals, because of the lattice peri- 
odicity, a roughening transition from a smooth to rough 
surface, occurs at some temperature. The nature of the 
growth of a crystal close to equilibrium expectedly de- 
pends on whether the surface is smooth or rough. One 
can also think of a crystal growth process which is far 
away from equilibrium by subjecting it to an external 
drive, for instance by random deposition of particles on 
the surface. The roughening that occurs in the nonequi- 
librium case is called kinetic roughening. Is the na- 
ture of the surface any different in kinetic roughening? 
Crystals are definitely not the only example of growth 
processes; some other examples of such nonequilibrium 
growths would be the growth of bacterial colonies in a 
petri dish, sedimentation of colloids in a drop, the forma- 
tion of clouds in the upper atmosphere, and so on. Note 
the large variation of length scales of these problems. In 
many such examples it is difficult if not impossible to 
think of an equilibrium counterpart. 



Scale invariance in interface fluctuations implies that 
fluctuations look statistically the same when viewed at 
different length scales. A quantitative measure of the 
height fluctuation (height measured from an arbitrary 
base) is provided by the correlation function 



C(x, t) = {[h{x + xo, < + to) - /i(xo, io)]^), 



(1) 



where x and t denote the d dimensional coordinate on the 
substrate and time respectively. The averaging in Eq. (|l|) 
is over all Xg, and, by definition, C(x, t) is independent 
of the choice of the arbitrary base. In simple language, 
scale invariance then means that when the system is, say, 
amplified by a scaling x ^ bx and t — > b^t, the height 
fluctuations reveal the same features as the original, upto 
an overall scale factor. Quantitatively, there exists a gen- 
eralized scaling 



C{x,t) = b-^^C{bx,bH), 



(2) 



where 6 is a scale factor, and x ^'Hd z are known as the 
roughening and dynamic exponents which are also uni- 
versal. As a direct consequence of (^), a scaling form for 
C{x, t) can be obtained by choosing b = 1/x 



C{x,t) 



C{t/x^), 



(3) 



a form that also explains the origin of the name "dynamic 
exponent" for z. The power law behaviour (as opposed to 
say exponential decay) of the correlation function implies 
absence of any scale, neither in space nor in time. All the 
underlying length scales required to deflne the problem 
dropped out of the leading behaviour in Eq. ^ Such a 
scale invariance is one of the most important features 
of equilibrium phase transitions and is observed when a 
parameter, say the temperature, approaches its critical 
value. However, here there is no special tuning param- 
eter; the scale invariance appears from the interplay of 
competing processes which in the simplest case can be 
the surface tension and noise present due to inherent ran- 
domness in the growth. There can be, of course, more 
complex events like a phase transition between surfaces 
with different roughness but scale invariance (not only of 
the correlation function but of any physical quantity) is 
generically preserved in all these surfaces. 
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It is worth emphasizing the enormous simpHfication 
that occurs in the scahng description. It is only a very 
few quantities that define the asymptotic behaviour of 
the system. Consequently, the idea of studying the uni- 
versal aspects of growth processes is to classify and char- 
acterize the various universality classes as determined by 
the exponents, e.g., x ^nd z, the scaling function and if 
necessary certain other important universal quantities. 

At this point, it might be helpful to compare the equi- 
librium and noncquilibrium cases again. In equilibrium, 
thanks to thermal energy (or "random kicks" from a heat 
reservoir), all configurations of a system are accessible 
and do occur, but no net flow of probability between 
any two states is expected (called "detailed balance"). 
Consequently, the knowledge of the states (and the en- 
ergies) of a system allows one to obtain the thermody- 
namic free energy by summing over the Boltzmann fac- 
tors exp(— iJ/fcsT), where E is the energy of the state, T 
the temperature and ks is the Boltzmann constant. In a 
noncquilibrium situation, either or both of the above two 
conditions may be violated, and the framework of pre- 
dicting the properties of a system from free energy is not 
necessarily available. A dynamic formulation is needed. 
By assigning a time dependent probability for the system 
to be in a configuration at a particular time, one may 
study the time evolution of the probability. The equilib- 
rium problem can be viewed from a dynamical point also. 
This description must give back the Boltzmann distribu- 
tion in the infinite time steady state limit. This is the 
Fokker Planck approach. The probabilistic description 
comes from the ensemble picture where identical copies 
of the same system exchange energy with the bath inde- 
pendently. An alternative approach which finds easy gen- 
eralization to the noncquilibrium cases is the Langevin 
approach where one describes the time evolution of the 
degrees of freedom, in our example /i(x, i), taking care of 
the random exchange of energy by a noise. The dynamics 
we would consider is dissipative so that the system in ab- 
sence of any noise would tend to a steady state. However 
for it to reach the equilibrium Boltzmann distribution in 
the presence of noise, it is clear that the noise must sat- 
isfy certain conditions (Einstein relation) connecting it 
to the system parameters. The nonequilibrium case does 
not have any thermodynamic free energy as a guiding 
light and therefore, there is no requirement to reach the 
Boltzmann distribution. In the Langevin approach, the 
noise term can be completely independent. In the equi- 
librium case the Langevin equation will be determined by 
the Hamiltonian or the free energy of the system, but for 
nonequilibrium cases there might be terms which can- 
not be obtained from Hamiltonians. Since for t oo, 
the probability distribution for equilibrium cases attains 
the Boltzmann distribution, the roughness exponent x 
is determined even in dynamics by the stationary state 
while the details of the dynamics is encoded in the dy- 
namic exponent z. In other words the two exponents x 



and z are independent quantities. In the nonequilibrium 
case, there is no compulsion to reach any predetermined 
stationary state and therefore the surface roughness is 
related to the growth, i.e., x and z need not be indepen- 
dent. We see below that there is in fact a specific relation 
connecting these two exponents. 

The existence of scale invariance and universal expo- 
nents implies that as far as the exponents are concerned 
the theory should be insensitive to the microscopic de- 
tails or in other words one may integrate out all the small 
length scale features. The universal exponents come out 
as an output of this process of coarse graining of say the 
Langevin equation, followed by a length rescaling that 
brings the system back to its original form. The new 
system will however have different values of the parame- 
ters and one can study the flow of these parameters in the 
long length and time scale limit. This is the basic idea 
behind the renormalization group (RG). In this approach 
the importance of an interaction or a term is judged not 
by its numerical value but rather by its relevance. One 
may start with any physically possible process in dynam- 
ics and see how it appears as the length scale or resolu- 
tion changes. For large length scales, one is left only with 
the relevant terms that grow with length, and marginal 
terms that do not change; the irrelevant terms or interac- 
tions that decay with length scale automatically drop out 
from the theory. The exponents are determined at the 
fixed points of the flows in the parameter space. These 
fixed points, which remain invariant under renormaliza- 
tion, characterize the macroscopic or asymptotic behav- 
ior of the system. Clearly from this view point of renor- 
malization group, one can explain why the microscopic 
details can be ignored and how the idea of "universality" 
emerges. All systems whose dynamical behaviour would 
flow to the same fixed point under RG transformation 
will have identical scaling behaviour. The various uni- 
versality classes can then be associated with the various 
fixed points of RG transformations, and phase transi- 
tions or criticality with unstable fixed points or special 
flows in the parameter space. An RG approach therefore 
seems rather natural and well suited for studying any 
scale invariant phenomena in general, growth problems 
in particular. Quite expectedly, the modern approach to 
growth problems is based on these views of RG. 

For a quantitative discussion, we consider two simple 
equations that, from the historical point of view, played 
a crucial role in the development of the subject in the 
last two decades. 

A simple Langevin equation describing the dynamics 
of a surface is the Edwards- Wilkinson (EW) equation ||] 

^ = „W^h + r^{x,t), (4) 

where x represents in general the coordinate on the d 
dimensional substrate, v is the coefficient of the diffu- 
sion term trying to smoothen the surface and 77 is the 
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Langevin noise which tries to roughen the surface 
One may add a constant current c to the right hand 
side, but by going over to a moving frame of reference 
{h ^ h + ct) one recovers Eq. ^. The noise here is 
chosen to have zero mean and short range correlation as 

(7?(x,t)77(x',i')) = 2£)(5(x-x')(5(t-t')- One of the impor- 
tant assumptions in this equation is that the surface is 
single- valued and there are no overhangs. One can solve 
exactly just by taking the Fourier transform and ob- 
tain the exponents x — (2 — (i)/2 and z = 2. That the 
dynamic exponent z is 2 follows from the simple fact that 
the equation involves the first derivative in time but a 
second derivative in space. The surface is logarithmically 
rough at d = 2. For d > 2, the fluctuations in the height 
are bounded and such a surface is more or less flat, better 
called "asymptotically flat" . From the growth equation 
one can also derive the stationary probability distribution 
for the height /i(x) which takes the form of a Boltzmann 
factor P(/i(x)) cx exp[— (i//iD) / (V/i)^(i'^a;] resembling an 
equilibrium system at a temperature given hy D — ksT. 
This is the Einstein relation noise should satisfy to re- 
cover equilibrium probability distribution. Conversely, 
given a hamiltonian of the form J {'Vh)^d'^x, the equi- 
librium dynamics will be given by Eq. ^ with D deter- 
mined by the temperature. Nevertheless, if we do not 
ascribe any thermal meaning to D, Eq. ^ is good enough 
to describe a nonequilibrium dynamic process as well. 
Such a nonequilibrium growth will have many similarities 
with equilibrium processes, differing only in the origin of 
the noise, e.g. the expected symmetry h — > —h with 
(h) = in equilibrium will be preserved in the nonequi- 
librium case also. The growing surface with a correlation 
C(x,t) =1 X |(2-rf)/2 C{t/\x\^) will be similar in both 
cases for d < 2. 

A genuine nonequilibrium process will involve breaking 
the up-down symmetry which in equilibrium foUws from 
detailed balance. It should therefore be represented by a 
term involving even powers of h. We already saw that a 
constant current (zeroth power) does not add anything 
new. Since the origin in space or time or the position 
of the basal plane should not matter, the first possible 
term is (V/i)^. By looking at the geometry of a rough 
surface, it is easy to see that such a term implies a lateral 
growth that would happen if a deposited particle sticks 
to the first particle it touches on the surface. One gets 
the Kardar-Parisi-Zhang (KPZ) equation |^ 

^^^T7^h+^{Vhr+7^{^,t). (5) 

As a consequence of its mapping to the noisy Burger's 
equation, to the statistical mechanics of directed poly- 
mer in a random medium and other equilibrium and 
nonequilibrium systems, the KPZ equation has become a 
model of quite widespread interest in statistical mechan- 
ics. Though we focus on growth problems in this paper, 
the KPZ equation is also applicable in erosion processes. 



Taking a cue from the development in understanding 
the growth phenomena through the KPZ equation, a vast 
class of simulational and analytical models have evolved 
to explain different experimentally observed growth pro- 
cesses. Diverse technical tools ranging from simula- 
tions with various dynamical rules to different versions 
of renormalization group techniques, mode coupling the- 
ory, transfer matrix techniques, scaling arguments have 
been employed to understand kinetic roughening. In this 
review we attempt to provide an overview of this phe- 
nomenon of roughening of a growing surface. It is almost 
beyond the scope of this review, to describe in detail var- 
ious models and their experimental relevance. Rather we 
focus our attention on a few examples which may broadly 
represent a few different routes along which research has 
continued. 

The plan of this article is as follows. In the next sec- 
tion we focus on the KPZ equation and its renormaliza- 
tion group description. We also point out the connection 
of the KPZ equation to some other problems of physics. 
In section III a more generalized growth mechanism in- 
volving nonlocal interactions is presented. Section IV is 
devoted to the progress in understanding the roughening 
and super roughening transitions which appear in a very 
distinct class of models involving lattice pinning poten- 
tial. 

II. KPZ EQUATION AND MORE 

Let us first look at the origin of the various terms in 
Eq. ^ and ^ In both the equations, the noise term rep- 
resents random deposition, the fluctuation around the 
steady value. As already mentioned a steady current can 
be removed from the equation by going over to a mov- 
ing frame. The term involving second derivative of h can 
represent either of two processes. It could be a surface 
tension controlled diffusion process, in which a particle 
comes to the surface and then does a random walk on the 
surface to settle at the minimum height position thereby 
smoothening the surface. An alternative interpretation 
would be that there is desorption from the surface and 
the process is proportional to the chemical potential gra- 
dient. The chemical potential of the particles on the sur- 
face cannot depend on h or gradient of h because of its 
independence of the arbitrary base or its tilt. The chem- 
ical potential then is related to the second derivative of 
h. This also has a geometric meaning that W^h is re- 
lated to the local curvature. The larger the curvature, 
the higher is the chance to desorb because of a lesser 
number of neighbours. In the KPZ equation, the non- 
linear term represents lateral growth. The diffusion-like 
term can then be thought of either (a) as an alternative 
that a particle coming to the surface instead of sticking 
to the first particle it touches, deposits on the surface 
and then diffuses, or (b) as a random deposition process 
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with desorption. In either case, the noise term tends 
to roughen the surface, the diffusion term, of whatever 
origin, smoothens it while the nonhnear term leads to a 
laterally growing surface. Even if the smoothening linear 
term is not present, renormalization group or the scaling 
argument indicates that such a term is generated on a 
large length scale. 

The KPZ equation has a special symmetry not present 
in the Edwards- Wilkinson case. This is the tilt symmetry 
(often called Galilean invariance - a misnomer, though, in 
this context). If we tilt the surface by a small angle, then 
with a reparametrization h' — h + e ■ x and x = x + Xet' 
and t = t' , the equation remains invariant for small e. 
This transformation depends only on A the coefhcient of 
the nonlinear term and fails for A = 0. Since this tilt sym- 
metry is to be maintained no matter at what lengthscale 
we look at, A must be a renormalization group invariant. 

Let us now perform a length rescaling analysis. Under 
a change of scale as a; ^ bx, t — > and h b^h, KPZ 
equation transforms as 

bx-^^ = „b^-^V'h + ^b^x-^{\/hf + b-'^^^-'/^rj, (6) 

where the noise correlation has been used to obtain the 
scaling of the noise term. Therefore under this scale 
transformation different parameters scale as — > 
D b'^-'^-'^^D and A b^+'--'^\. For A = 0, the equa- 
tion remains invariant provided z — 2 and % = (2 — d)/2. 
These are just the exponents one expects from the EW 
model. (Such surfaces with anisotropic scaling in dif- 
ferent directions like x and h are called self-afhne. ) 
Though we cannot predict the exponents from Eq. ^ 
when A 7^ 0, it does tell us that a small nonlinearity 
added to the EW equation scales with a scaling dimen- 
sion X + ~ 2. This term is always relevant in one di- 
mension, because it scales like 6^/^. This type of scal- 
ing argument also shows that no other integral powers 
of derivatives of h need be considered in Eq. |^ as they 
are all irrelevant, except (dh/dx)^ at d = 1, which, how- 
ever, detailed analysis shows to be marginally irrelevant. 
Based on this analysis, we reach an important conclusion 
that the nonequilibrium behaviour in one dimension, and 
in fact for any dimensions below two, would be distinctly 
different from the equilibrium behaviour. For dimensions 
greater than two, Edwards Wilkinson or equilibrium sur- 
faces, as already mentioned, are asymptotically flat with 
X = 0, 2 = 2, and so, a small nonlinearity is irrelevant 
because it will decay with b. In other words, the growth 
in higher dimensions for small A would be very similar 
to equilibrium problems because the EW model is stable 
with respect to a small perturbation with nonlinearity. 
The simple scaling argument does not tell us if the na- 
ture of the surface changes for large A for c? > 2, but 
an RG analysis shows that it does change. That the 
nonequilibrium growth is always different in low dimen- 
sions and in higher dimensions (greater than two), and 



that there will be a dynamic phase transition from an 
equilibrium-like to a genuine nonequilibrium behaviour, 
explains the source of excitement in this minimal KPZ 
equation, in the last two decades. 

If the nonlinear parameter A is to remain an invariant, 
i.e. independent of b in Eq. (^, then x-|-z = 2, a relation 
which need not be satisfied by the equilibrium growth. It 
is this relation connecting the two exponents of the scal- 
ing function of Eq. |^ that distinguishes nonequilibrium 
growth from equilibrium, the former requiring one less 
exponent than the latter one. We wonder if such an ex- 
ponent relation is generally true for all nonequilibrium 
systems. 

Though we are far away from a complete understand- 
ing of all the nuances and details of the KPZ equation, 
the renormalization group analysis has been very suc- 
cessful in identifying different phases, nature of phase 
transitions, and, in certain cases, relevant exponents. In 
brief, the various results obtained from renormalization 
group analysis are as follows. In one dimension and for 
d < 2, even a small nonlinearity, as already mentioned, 
being relevant in the RG sense, leads to new values of 
roughening and dynamic exponents, and is characterized 
by a RG fixed point. Beyond d = 2, there is a phase 
transition demarcating two different types of surfaces. A 
small nonlinearity is irrelevant around EW model and 
the surface is almost flat with x = and z = 2. A strong 
nonlinear growth, however, drives the system to a dif- 
ferent phase with rougher surface where x 7^ 0. Several 
aspects of this phase transition can be studied from RG 
but the strong A regime is still out of reach, because of 
the absence of any RG fixed point. 

The KPZ equation in d = 1 has distinct nonequilib- 
rium behaviour, and the scaling behaviour is the same no 
matter how small or large A is. More peculiar is the exis- 
tence of a stationary probability distribution of the height 
in one dimension which is the same as for the linear EW 
model. This is not just an accident but a consequence of 
certain subtle relations valid only in one dimension. We 
do not go into those issues here. The same stationary dis- 
tribution implies that the nonlinearity does not affect the 
stationary state solution, and x = 1/2. The two mod- 
els however differ in the dynamic exponent which, in the 
case of KPZ growth, has to satisfy x + ^ = 2. This leads 
to exact answer z — 3/2. Its significance can be grasped 
if we compare various known cases. For ballistic motion, 
distance goes linearly with time so that the dynamic ex- 
ponent is z = 1 while for diffusive motion or in quantum 
mechanics (e.g. a nonrelativistic free quantum particle), 
z = 2 as also the case for EW. Here is an example where 
the nonequilibrium nature of the problem leads to a com- 
pletely new exponent connecting the scaling of space and 
that of time. 
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A. Dynamic renormalization group analysis 

A dynamic renormalization group analysis is a more 
general approach applicable for dynamics which e.g. may 
be governed by the Langevin equation for the appropri- 
ate dynamical variable. For our problem it is easier to 
work in Fourier coordinates q, and ui conjugate to space 
and time. Long distance, long time implies q — > and 
w — > 0, and q can be taken as the inverse wavelength at 
which the height variable is probed. The magnitude of 
wave vector q varies from to A where the upper cutoff 
is determined by the underlying microscopic length scale 
like lattice spacing or size of particles etc. In the Fourier 
space, different Fourier modes in the linear EW model 
gets decoupled so that each /i(q, lj) for each (q, w) be- 
haves independently. It is this decoupling that allows the 
simple rescaling analysis of Eq. ^ or dimensional anal- 
ysis to give the correct exponents. For the KPZ equa- 
tion the nonlinear term couples heights of various wave- 
lengths and therefore any attempt to integrate out the 
large (q, tf) modes will affect h with low values of (q, w). 
This mixing is taken into account in the RG analysis 
which is implemented in a perturbative way. One thinks 
of the noise and the nonlinear term as disturbances af- 
fecting the EW-like surface. If we know the response of 
such a surface to a localized disturbance we may recover 
the full response by summing over the disturbances at 
all the points and times. However this disturbance from 
the nonlinear term itself depends on the height, requir- 
ing an iterative approach that generates successively a 
series of terms. By averaging over the noise, one then 
can compute any physical quantity. At this stage only 
degrees of freedom with q in a small shell e~'A < g < A 
is integrated out. In real space this corresponds to inte- 
grating out the small scale fluctuation. The contribution 
from this integration over the shell is absorbed by re- 
defining the various parameters i/, A and D. These are 
the coupling constants for a similar equation as but 
with a smaller cutoff Ae"'. A subsequent rescaling then 
restores the original cutoff to A. Following this proce- 
dure, the flow equations for different parameters v, D, 
and A can be obtained Using the exponent identity 
predicted from the Galilean invariance and the renormal- 
ization group invariance of v, the flow equations for all 
the parameters can be combined into a single flow equa- 
tion for A = X'^D/v^ (with A = 1). This is the only 
dimensionless parameter that can be constructed from 
A, V, D, and A, and it is always easier to work with 
dimensionless quantities. Its recursion relation is 

dX 2~d- 2d - 3-3 

where Kd is the surface area of a d-dimensional sphere 
divided by (27r)'^. The invariance of i/ under RG trans- 
formation implies z = 2 — KdX ^rj^, and the Galilean 



invariance provides the value of x = 2 — z once the value 
of z is known. To be noted here is that the dynamic ex- 
ponent is different from 2 by a term that depends on A 
coming from the renormalization effects. 

A few very important features are apparent from (|^). 
From the fixed point requirement dX/dl = 0, we find that 

2 

at = 1, there is a stable fixed point A = 2/Ki. At this 
fixed point z = 3/2 and % = 1/2 supporting the results 
predicted from the symmetry analysis. At d = 2, the 
coupling is marginally relevant indicating a strong cou- 
pling phase not accessible in a perturbation scheme. At 
d > 2, the flow equation indicates two distinct regimes, 
namely a weak coupling regime where A asymptotically 
vanishes leading to a flat EW phase with x = 0, z = 2, 
and a strong coupling rough phase, the fixed point of 
which cannot be reached by perturbation analysis. 

Owing to this limitation of the renormalization group 
analysis based on the perturbation expansion, the scaling 
exponents in this strong coupling phase cannot be deter- 
mined by this RG scheme. Different numerical methods 
yield z — 1.6 ai d — 2. The phase transition governed by 
the unstable fixed point of A is well under control with 
z = 2 ioT all d > 2. To explore the strong coupling phase, 
recently techniques like self-consistent mode coupling ap- 
proach, functional renormalization group etc have been 
employed, but even a basic question whether there is an 
upper critical dimension at which z will again become 2 
remains controversial. 



B. Relation with other systems 

The relation of the KPZ equation with other quite un- 
related topics in equilibrium and nonequilibrium statisti- 
cal mechanics is impressive. Here we provide a very brief 
account of these systems. 

Noisy Burgers equation: By defining a new variable 

V = V/i, we obtain an equation 

^ZJV^v + Av- Vw + f(x,i), (8) 

ot 

where the noise term / = V?7. The above equation 
represents the noisy Burgers equation for vortex free ( 

V X V = 0) fluid flow with a random force. This equa- 
tion is very important in studies of turbulence. The tilt 
invariance of the KPZ equation turns out to be the con- 
ventional Galilean invariance for the Burgers equation 
(for A = 1), and that's how the name stayed on. 

Directed polymer in a random medium: A directed 
polymer, very frequently encountered in different prob- 
lems in statistical mechanics, is a string like object which 
has a preferred longitudinal direction along which it is 
oriented, with fluctuations in the transverse direction. 
The flux lines in type II or high Tc superconductors 
are examples of such directed polymers in 3 dimensions, 
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while the steps on a vicinal or miscut crystal surface 
or the domain walls in a uniaxial two dimensional sys- 
tem are examples in two dimensions. The formal math- 
ematical mapping to such objects follows from a simple 
(Cole-Hopf) transformation of the KPZ equation using 
W{x,t) — exp[jj^h]. The Cole-Hopf transformation lin- 
earizes the nonlinear KPZ equation and the resulting lin- 
ear diffusion equation (or imaginary time Schroedinger 
equation) is identical to that satisfied by the partition 
function of a directed polymer in a random potential. For 
such random problems, one is generally interested in the 
averages of thermodynamic quantities like the free energy 
and we see that the noise averaged height (/i(x, t)) gives 
the average free energy of a directed polymer of length 
t with one end at origin and the other end at x. This 
is a unique example of a system where the effect of such 
quenched averaging of free energy can be studied without 
invoking any tricks (like the replica method). This has 
led to many important results and enriched our under- 
standing of equilibrium statistical mechanics. Recently, 
this formulation has been extended to study details of the 
properties of the random system near the phase transi- 
tion point and overlaps in lower dimensions . It turns 
out that one needs an infinite number of exponents to 
describe the statistical behaviour of the configurations of 
the polymer in the random medium We do not go 
into this issue as this is beyond the scope of this article. 

An interesting connection between the 1-1-1 dimen- 
sional KPZ equation and the equilibrium statistical me- 
chanics of a two dimensional smectic -A liquid crystals 
has been recently established by Golubovich and Wang 
. This relationship further provides exact approach to 
study the anomalous elasticity of smectic-A liquid crys- 
tals. 

Apart from these, there are a number of other relations 
between KPZ equation and kinetics of annihilation pro- 
cesses with driven diffusion, the sine-Gordon chain, the 
driven diffusion equation and so on. 

C. Beyond KPZ 

1. Conservation condition 

The situation encountered in Molecular-beam epitaxy 
(MBE) for growth of thin films is quite different than the 
mechanism prescribed by the KPZ equation In MBE 
the surface diffusion takes place according to the chemical 
potential gradient on the surface, respecting the conser- 
vation of particles. If the particle concentration does not 
vary during growth, then a mass conservation leads to a 
volume conservation and the film thickness is governed 
by an underlying continuity equation 

dh „ . 

7^ + V-j=ry, (9) 



where j is the surface diffusion current which states that 
the change of height at one point is due to flow into or out 
from that point. The current is then determined by the 
gradient of the chemical potential, and since the chemical 
potential has already been argued to be proportional to 
the curvature V^/i, the growth equation thus becomes a 
simple linear equation involving V^h which, like the EW 
model, is exactly solvable. Taking into account the effect 
of nonlinearity the full equation can be written as 

^ = -\7'WV'h+^iVhf]+r^ix,t), (10) 

where the noise correlation is {ri{x,t)r]{x' ,t')) = 
2DV^S{x — x')6{t — t'), if the noise also maintains con- 
servation (if it originates from the stochasticity of diffu- 
sion) or would be the same white noise as in the KPZ 
equation if the noise is from random deposition. It goes 
without saying that the exponents are different from the 
EW model even for the linear theory. The invariance of A 
in this case leads to a different relation between x and z. 
At the dimension of physical interest d — 2, this growth 
equation leads to an enhanced roughness than the KPZ 
case and may explain the results of experiments of high 
temperature MBE. 



2. Quenched noise 

A different type of generalization of the KPZ equation 
was to explore the motion of domain walls or interfaces 
in a random medium. In this case, the noise is not ex- 
plicitly dependent on time but on the spatial position 
and the height variable. Such a noise has been called 
quenched noise because the noise is predetermined and 
the interface or the surface moves in this random sys- 
tem. The simple features of the KPZ equation and the 
EW model are lost. Functional rcnormalization group 
analysis and numerical studies attempted to clarify the 
question of universality classes and details of dynamics in 
such cases. The important concept that emerged in this 
context is the depinning transition so that the surface re- 
mains pinned by the randomness until the drive exceeds 
a certain critical value. Interface depinning is an exam- 
ple of a nonequilibrium phase transition. The velocity of 
the surface near this depinning transition also has crit- 
ical like behaviour with long range correlations. Below 
the threshold, the dynamics is sluggish, while just above 
the threshold , the velocity is in general not proportional 
to the drive but obeys a power law with a universal ex- 
ponent. For a very strong drive (or large velocity of the 
interface) the moving surface encounters each site only 
once, and therefore the noise is effectively like a space- 
time dependent noise rather than the quenched one. The 
nature of the surface would then be like KPZ. 
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3. Coloured noise 

In the previous section we discussed the KPZ equation 
with white noise. If the noise is coloured in the sense 
that there is correlation in space or time or both, the 
universal behaviour, the phase transitions and the prop- 
erties are different but still can be studied by the same 
RG technique. Several aspects of the problem especially 
the role of noise correlation have been explored | |lO[ |. 

All of the above seem to suggest that if there is no 
conservation law, then the KPZ equation is the equation 
to describe any nonlinear or nonequilibrium growth pro- 
cess and all phenomena can be put in one of the known 
universality classes. However, experimentally KPZ expo- 
nents seem to elude us so far Since results are 
known exactly in one dimension, special one dimensional 
experiments were conducted like paper burning, inter- 
face motion in paper, colloid suspension etc, but KPZ 
exponents have not been seen. In the colloid experiment 
PI the surface formed by the depositing colloids on the 
contact line (d = 1) between the colloid latex film and a 
glass slide was measured from video images. This method 
yields x = -71 but cannot determine the dynamic expo- 
nent. A recent analysis of tropical cumulus cloud in 
the upper atmosphere, from satellite and space shuttle 
data from 0.1 to 1000 Km, seems to agree with the KPZ 
results in d = 2. 



III. KINETIC ROUGHENING WITH 
NONLOCALITY 



In spite of a tremendous conceptual and quantitative 
success of KPZ equation in describing the nonequilib- 
rium growth mechanism, the agreement with experimen- 
tally observed exponents is rather unsatisfactory. One 
wonders whether there is any relevant perturbation that 
drives the systems away from the KPZ strong coupling 
perturbation. One goal of this section is to point out that 
indeed there can be longrange interactions that may give 
rise to non-KPZ fixed points. 

Many recently studied systems involving proteins, col- 
loids, or latex particles the medium induced interactions 
are found to play an important role [T^ . This nonlo- 
cal interaction can be introduced by making a modifica- 
tion of the nonlinear term in the KPZ equation. Taking 
the gradient term as the measure of the local density 
of deposited particles, the long range effect is incorpo- 
rated by coupling these gradients at two different points. 
The resulting growth equation is a KPZ-like equation 
with the nonlinear term modified as ^ / (ir'V(r') 
V/i(r -I- r', t) • V/i(r — r',i). For generality we take V(r') 
to have both short and long range parts with a specific 
form in Fourier space as V(A:) = Aq -t- Xpk~'' such that in 



to observe whether the macroscopic properties are gov- 
erned by only Aq and hence KPZ like or the behavior is 
completely different from KPZ due to the relevance of Ap 
around the KPZ fixed points. 

A scaling analysis as done in Eq. (||) clearly indicates 
different scaling regimes and the relevance of Aq and Ap 
for (i < 2 at the EW fixed point. For any \p{^ 0) with 
p > 0, the local KPZ theory (i.e. Ap = and x + 2 = 2) is 
unstable under renormalization and a non-KPZ behavior 
is expected. For 2 < d < 2 -f 2p, only Ap is relevant at the 
EW fixed point. The exponents of the non-KPZ phases 
can be obtained by performing a dynamic renormaliza- 
tion group calculation By identifying the phases 
with the stable fixed points, we then see the emergence 
of a new fixed point where the long range features dom- 
inate (x + ^ = 2 -f p). Most importantly, at c? = 2, the 
marginal relevance of A is lost and there is a stable fixed 
point (LR) for p > .0194. 

On the experimental side, there are experiments on col- 
loids with X = 0.71 which is the value also obtained from 
paper burning exponents. For colloids, hydrodynamic in- 
teractions are important. Similar longrange interactions 
could also play a role in paper burning experiment due to 
the microstructure of the paper. With this x our expo- 
nents suggest p = —0.12 at d = 1 at the long range fixed 
point. Further experiments on deposition of latex parti- 
cles or proteins yielding the roughness of growing surface 
have not been performed. Probably such experiments 
may reveal more insights on this growth mechanism. 

More recently, the effect of coloured noise in presence 
of nonlocality has been studied ||l^ and the nature of 
the phases and the various phase transitions clarified. A 
conserved version of the nonlocal equation has also been 



considered and it shows rich behaviour |16 



the limit Xo 



KPZ results are retrieved. The aim is 



IV. ROUGHENING TRANSITION IN 
NONEQUILIBRIUM 

It is interesting to study the impact of equilibrium 
phase transitions on the nonequilibrium growth of a sur- 
face. This is the situation observed experimentally in 
growth of solid '^He in contact with the superfluid phase 
|]l7| . There is an equilibrium roughening transition at 
Tfi — 1.28K. For T > Tr the growth velocity is hnear in 
the driving force F (chemical potential difference), but 
for T < the velocity is exponentially small in the in- 
verse of the driving force. For infinitesimal drive , the 
mobility which is the ratio of the growth velocity and F 
vanishes with a jump from a finite value at the transi- 
tion. With a finite force the transition is blurred and the 
flat phase below in equilibrium becomes rougher over 
large length scale. 

The equilibrium roughening transition is an effect of 
discrete translational symmetry of the lattice. The equi- 
librium dynamics in this case is essentially governed by 
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the Langevin equation 

^^KS/^h{r,t)-Vsm[^ h{r,t)] + ar,t), (11) 

where the sin term favours a periodic structure of spao- 
ing a. Extensive investigations have been done on this 
equihbrium model. At low temperature this periodic po- 
tential is relevant and it ensures that minimum energy 
configuration is achieved when is an integer multiple 
of lattice periodicity. In this phase the surface is smooth 
and the roughness is independent of length. In the high 
temperature phase the equilibrium surface is thermally 
rough and the roughness is logarithmic 



C{L,T)^\n[Lf{r/L^-)]. 



(12) 



The critical point is rather complicated and goes by the 
name of Kosterlitz-Thouless transition, first discussed in 
the context of defect mediated transitions in two dimen- 
sional XY magnets . 

For a nonequilibrium crystal growth problem, one 
needs to introduce the KPZ nonlinear term in (11). 



There is no longer any roughening transition. The fact 
that away from equilibrium the roughening transition is 
blurred is manifested by the domination of the nonlinear 
term and the suppression of the pinning potential in the 
asymptotic regime ||l^ . 

A very nontrivial situation arises when the surface con- 
tains quenched disorder which shifts the position of the 
minima of the pinning potential in an arbitrary random 
fashion |20| , pT[ |. In this case there is a new phase transi- 
tion which is drastically different from the equilibrium 
roughening transition. This transition is called super 
roughening. Above the transition temperature i.e. for 
T > Tsr-, the surface is logarithmically rough as it is in 
the high temperature phase of the pure problem. How- 
ever in the low temperature phase i.e. for T < Tsr, the 
surface is no longer flat and is even rougher than the 
high temperature phase. Recent numerical treatments 
suggest that the surface roughness behaves as (InL)^. 
In the nonequilibrium situation the linear response mo- 
bility vanishes continuously at the transition tempera- 
ture unlike the jump discontinuity in the pure case. A 
general treatment with a correlated disorder elucidates 
the connection between the roughening and super rough- 
ening transition and one observes that the roughening 
turns into a super roughening transition if the disorder 
correlation decays sufhciently fast. Away from equilib- 
rium, the super roughening transition is essentially dom- 
inated by the KPZ nonlinearity and instead of the loga- 
rithmic roughness, an asymptotic power law behavior of 
the roughness is found over all temperature range. 

In a similar situation in the nonequilibrium case, one 
needs to study the role of the KPZ nonlinearity with 
long range disorder correlation [ p2[ . A functional renor- 
malization scheme with an arbitrary form of the disorder 



correlation turns out to be useful, though a detailed so- 
lution is not available. It is found that the flow of the 
KPZ nonlinearity under renormalization, with power law 
form of the disorder correlation, is such that it decays 
with length. This implies that nonequilibrium feature 
does not set in over a certain length scale. Over this 
scale one would then expect usual roughening transition. 
However there is a generation of a driving force due to 
the nonlinearity, and the growth of this force with length 
scale would invalidate use of perturbative analysis. For 
large length scales, one expects a KPZ type power law 
roughness of the surface. Nevertheless, the initial decay 
of the nonlinearity with the length scale due to the long 
range correlation of the disorder is an interesting conclu- 
sion that seems to be experimentally detectable. 



V. REMARKS 

In this brief overview, we attempted to focus on the dif- 
ference between equilibrium and nonequilibrium growth 
problems with an emphasis on the scaling behaviour and 
renormalization group approach. Many details with ref- 
erences to pre-1995 papers can be found in Ref. [|l|,D, 
which should be consulted for more detailed analysis. 
Though the success story of the KPZ equation is rather 
impressive, there are still many unresolved, controversial 
issues. In fact for higher dimensions, the behaviour is not 
known with as much confidence as for lower dimensions. 
Developments in this direction are awaited. 

Note added in Proof 

(1) The growth mechanism of metal-organic films de- 
posited by the Langmuir-Blodgett technique has been 
studied in Ref. 23 by X-ray scattering and atomic force 
microscopy. The results have been interpreted by a com- 
bination of 1-dimensional EW equation (Eq. 4 ) and 
2-dimensional linear conserved equation (Eq. 10) with 
conserved noise. 

(2) For effects of nonlocality in equilibrium critical dy- 
namics, see Ref. 24. 
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